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ABSTRACT 


An  approximate  procedure  is  used  to  predict  the  effect  on  air  ioni¬ 
zation  of  small  amounts  of  cesium  and  sodium  in  the  ablating  surface 
material  of  pointed  cones.  Ionization  is  assumed  to  occur  due  to  finite  rate 
gas  phase  chemical  reactions  in  the  laminar  boundary  layer.  An  11  species 
(N^.  O^,  N,  O,  NO,  Cs  or  Na,  NO^,  Cs^  or  Na^,  ,  O  ,  e  )  16  reaction 
chemical  system  is  assumed. 
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INTRODUCTION 

For  slender  reentry  bodies  in  hypersonic  flight,  the  highest  tem¬ 
peratures  in  the  flow  field  may  be  attained  in  the  boundary  layer.  Pre¬ 
diction  of  air  ionization  must  then  be  based  on  solution  cf  the  boundary 
layer  equations.  Ionization  in  air  is  a  complex  process  involving  a  large 
number  of  species  and  chemical  reactions.  For  this  reason,  numerical 
computation  procedures  are  necessary  to  make  such  predictions,  particularly 

if  the  effect  of  finite  chemical  rates  must  be  taken  into  consideration. 

1,2  3 

Blottner  and  others  have  made  such  computations  of  ionization  in  hyper¬ 
sonic  laminar  boundary  layers  for  pure  air.  The  question  arises  whether 
some  foreign  species  that  contaminate  the  air  in  the  boundary  layer  from 
an  ablating  or  evaporating  solid  surface  may  not  retard  or  promote  the 
ionization  process.  Easily  ionizable  gases,  such  as  alkali  metal  vapors, 
suggest  themselves  as  “seeding"  materials  to  promote  the  ionization  process 
for  various  applications.  The  purpose  of  the  present  investigation  is  to 
predict  the  effects  on  air  ionization  of  small  amounts  of  cesium  or  sodium- 

added  to  an  ablating  inert  material  that  enters  the  boundary  layer  on  pointed 

4 

cones.  As  in  other  studies  of  boundary  layers  with  complicated  chemistry  , 
the  effect  of  trace  species  on  gross  thermodynamic  and  fluid  mechanic 
properties  of  the  boundary  layer  is  not  considered.  For  simplicity  it  was 
assumed  that  the  ablating  inert  material  has  the  properties  of  air,  thus 
neglecting  any  complicating  effects  of  hydrocarbons  on  the  ionization  process. 
ANALYSIS 

The  laminar  boundary  layer  equations  on  pointed  cones  are  best 
written  in  transformed  V  ,  ^  coordinates^' 
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If  velocity  u  is  expressed  as 
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and  the  continuity  equation  is  used  to  express  v  in  terms  of  f  the  momentum 
and  energy  equations  taike  on  the  following  form: 
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The  species  conservation  equation  is: 
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where  the  flux  velocities  q,  are  related  to  the  gradients  by  equations  of  the 
type: 
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Plasma  neutrality  will  be  assumed  due  to  the  small  Debye  shielding 
distance  at  the  conditions  that  are  of  practical  interest.  Thus,  the  diffusion 
of  ions  will  be  of  the  ambi-polar  type  and  free  diffusion  of  electrons  need 
not  be  considered.  However,  the  usual  expressions  for  ambipolar  diffusion 
must  be  modified  due  to  the  simultaneous  presence  of  electrons  and  both 
positive  and  negative  ions.  For  charged  particles  equation  (7)  is  modified 
by  the  presence  of  a  body  force  caused  by  a  microscopic  electric  field  due 


to  local  charge  separation: 
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Equation  (8)  is  valid  for  all  charged  particles  including  electrons.  For  the 
assumed  neutral  plasma,  the  electron  and  ion  concentrations  and  fluxes  are 
related  by  charg-  concentration  laws: 
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Equation  (9)  makes  the  electron  diffusion  equation  (8)  superfluous;  the 
latter  can  then  be  used  to  elir^inate  electric  field  E  from  the  remaining 
equations  (8).  The  result  is: 
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where  the  total  number  of  species  n  in  (10)  and  subsequent  equations  does 
not  include  electrons.  Due  to  the  large  mobility  of  electrons,  the  ratio 
&..  /&  .  will  be  a  small  number; 
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Eliminating  te^ms  in  (10)  proportional  to  the  small  ratio  in  (11)  one  obtains 
the  generalized  ambipolar  diffusion  law  for  positive  ions  ; 
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For  negative  ions,  the  simplification  of  (10)  is  not  so  straightforward  be¬ 
cause  the  sum  on  the  right  hand  side  and  the  first  sum  on  the  left  hand  side  will 
both  be  small,  so  that  the  small  term  involving  q^  cannot  any  more  be 
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considered  insignificant  compared  to  the  remaining  terms.  Thus,  for 
negative  ions. 
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A  corollary  of  the  foregoing  development  is  tha.  is  significantly  larger 
than  q.  .  This  fact  can  be  utilized  in  conservation  law  (9),  and  the  result 
substituted  into  (13); 
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Equation  (10)  represents  the  most  general  form  of  ambipolar  diffusion;  it 
is  based  on  the  assumption  of  no  macroscopic  charge  separation.  Equations 
(12)  and  (14)  are  based  on  three  additional  assumptions: 

(a)  &..<<&  .  (eq.  1 1) 

ij  €J  ^ 

(b)  ,  =5;  1.0  for  all  species  except  electrons 

ij  kl 

(c)  (c.  m  /c  m.)  constant 

1  c  C  1 

There  assumptions  are  reasonable  for  chemi-ionization  in  the  boundary 
layer  where  all  species  in  the  mixture  have  similar  molecular  weights, 

K  the  further  simplification  of  a  "constant  Lewis  number"  is  made. 


1.  e. : 


&  =  &  ,  &  .  =  & 
ij  f  J  ^ 


equations  (7),  (12)  and  (14)  reduce  to; 
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Results  (16),  together  with  (9),  represent  the  generalizations  of  Pick's  law 
in  the  boundary  layer  for  ambipolar  diffusion  in  a  many-ion  mixture. 

Thermodynamic  and  chemical  kinetic  relationships  are  required  to 
complete  the  system  of  equations: 


=  constant 


w.  =  w.  (p,  T,  c.)  (17) 

Except  for  the  role  of  negative  ions,  the  foregoing  system  of  equations 
is  derived  and  discussed  in  detail  in  reference  5;  numerical  results  for  the 
pure  air  system  are  given  in  references  1,  Z  and  5;  details  of  thermodynamic 
transport  and  chemical  properties  that  were  employed  are  given  in  reference 
5,  It  is  noteworthy  that  for  a  pointed  cone  there  was  no  substantial  difference 
between  numerical  results  of  reacting  and  ionizing  air  boundary  layer  ob¬ 
tained  for  the  general  non-similar  case  by  un  exact  implicit  finite  difference 
scheme  \  or  the  so-called  locally  similar  solution^,  which  neglects  all 
derivatives  with  respect  to 
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The  physico-chemical  processes  that  take  place  in  the  laminar 
hypersonic  boundary  layer  on  cones  can  be  ascertained  both  from  the  species 
conservation  equation  (6)  and  the  numerical  results  obtained  for  pure  air. 

The  quantity  ^  (du^/dx)  ^  that  determines  how  large  the  chemical  mass 
generation  rate  terms  are  in  the  differential  equation,  increases  from  0  at 
the  tip  of  the  cone  to  large  values  far  downstream.  Correspondingly,  the 
numerical  results  in  references  1,  2,  and  5  indicate  that  near  the  tip  the 
boundary  layer  composition  is  essentially  frozen,  and  as  distance  increases, 
and  the  boundary  layer  thickness  grows,  a  gradual  relaxation  towards 
equilibrium  composition  in  the  boundary  layer  takes  place.  The  principal 
chemical  activity  occurs  in  a  narrow  high  temperature  region  of  the  boundary 
layer  from  where  the  reaction  rc-o(f  cts  diffuse  towards  the  catalytic  wall 
surface  and  towards  the  free  st.  f  or  a  reasonably  cold  wall  surface 

and  for  slender  cones  where  the  inviscid  flow  temperature  is  low,  the  con¬ 
centrations  of  the  reaction  products  at  a  catalytic  wall  and  in  the  free  stream 
will  both  be  insignificant.  This  makes  it  possible  to  introduce  the  assumption 
of  "similar"  composition  profiles  for  all  reaction  products  in  the  boundary 
layer  by  means  of  a  function  of  f]  and  ^  ,  i,  e,  : 

c.  (n.  f)  =  c,  F  (n.  C  )  (18) 

max 

for  all  c.  that  are  rea  ;tion  products.  Function  F  has  the  desired  properties 
of  vanishing  at  the  wall  and  far  from  the  wall  and  rising  to  a  sharp  peak  in 
the  high  temperature  region.  Assumption  (18)  is  the  basis  of  the  present 
investigation  of  estimating  the  effect  of  small  amounts  of  cesium  or  sodium 
in  the  ablating  inert  material  of  the  cone  surface. 

Let  three  assumptions  now  be  made  to  simplify  the  description  of 
the  inter -diffusion  of  chemical  species  in  the  boundary  layer: 

(a)  All  binary  diffusion  coefficients  are  equal  at  a  given  temperature 

(b)  For  the  purposes  of  diffusion  laws  all  species  have  the  same 
molecular  weights. 

(c)  The  resulting  Schmidt  number  is  a  constant  regardless  of 
temperature. 
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Subject  to  the  above  assumptions,  Pick's  law  is  recovered,  and  the 
species  conservation  equations  (6)  and  diflusion  equations  (6)  for  neutral 
species  can  be  combined: 
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The  total  mass  fraction  of  constituent  elements  in  the  boundary 
layer  can  be  related  to  the  mass  fraction  of  the  individual  species: 
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where  a .  is  the  number  of  atoms  of  elements  j  in  species  i,  and  the 
m  -  s  are  the  respective  molecular  (or  atomic)  weights.  Since  elements 
are  conserved  in  a  chemical  reaction,  i.  e.: 
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an  equation  for  conservation  of  elements  may  be  written  by  neglecting  trace 
amounts  of  ionized  species  and  by  applying  summation  (20)  to  equations  (19): 
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For  pointed  cones  /3  =  O,  so  that  comparison  of  (3)  and  (22)  reveals 
that  a  "Reynolds  analogy"  can  be  established  for  the  element  concentrations 
by  a  simple  transformation  wh-'-h  incorporates  the  Schmidt  number  into  the 
coordinate  Tj  ; 
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provided  elemental  concentration  c'^  (0,  ^)  is  independent  of  ^  or  "local" 
similarity  is  assumed.  (It  will  be  shown  below  that  for  realistic  injection 
rate  distribution  the  latter  assumption  is  necessary.)  What  the  concentrations 
of  various  elements  are  at  the  wall  remains  to  be  determined  from  the 
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boundary  conditions.  The  general  boundary  conditions  for  species  diffusing 
to  and  from  a  catalytic  wall  where  material  is  injected  are: 

r  .  .  n  m.  .  j 

V  (0,x)  =  E  - ^  a"!  c.  q.!  (24) 
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where  v  is  an  injection  velocity  corresponding  to  the  mass  injecti'^n  rate 
w 

and  density  on  the  wall  surface;  and  c'^  the  mass  fraction  of  the  j-th  element 

w 

in  the  injected  gas.  Boundary  condition  (24)  accounts  for  the  fact  that  the 
total  amount  of  element  j  in  the  species  diffusing  to  and  from  the  wall  changes 
due  to  the  enrichment  by  the  gas  that  is  ejected  from  the  surface  (c"^  )  as  it 
mixes  with  the  boundary  layer  fluid  (c"^  (0,  x).  Using  the  transformed 
coordinate  system  and  the  approximations  implicit  in  (16),  boundary 
condition  (24)  can  be  rewritten: 
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Using  (23)  in  conjunctions  with  (2  5)  one  may  solve  for  the  elemental  com¬ 
position  at  the  wall: 
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where  I  is  the  injection  parameter; 
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and  the  second  form  of  I  in  (27)  is  based  on  constant  wall  and  free  stream 

properties  for  a  cone.  It  can  be  seen  from  (26)  and  (27)  that,  except  for 

special  conditions,  (0,  ^)  will  be  a  function  of  ^  .  Expression  (23)  must 

then  be  justified  by  the  assumption  of  local  similarity.  As  mentioned  before, 

5 

for  pure  air  this  assumption  was  found  reasonably  valid  ,  and,  unless  there 
are  sudden  variations  of  c^'^  and  1  along  the  body,  it  still  appears  reasonable 
for  contaminated  air. 
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Now  finally  the  problem  of  injecting  an  inert  gas  containing  small 
amounts  of  cesium  or  sodium  can  be  considered.  Let  the  injected  inert  gas 
have  the  properties  of  air,  and  let  the  amount  of  cesium  or  sodium  be  suf¬ 
ficiently  small  so  that  the  gross  thermal  and  fluid  mechanic  properties  of 

7 

the  flow  remain  unchanged  by  the  addition.  The  air  chemistry  used  in 
reference  1,  2  and  5  can  then  be  used  to  calculate  properties  of  the  boundary 
layer  with  air  injection  for  a  given  flight  condition.  This  system  considers 
NO  as  the  sole  source  of  ionized  particles.  Let  the  total  cesium  or  sodium 
concentration  be  determined  by  (26)  and  (23),  Furthermore,  let  the  additional 
ionization  products  due  to  the  presence  of  cesium  be  related  to  the  NO  ions 
by  the  "similarity  law"  of  equation  (18).  Differential  equation  (19)  can  now 
be  considered  for  the  ionization  products  at  the  point  where  the  ion  con¬ 
centrations  are  the  maximum,  (i,  e.  SF/Bt)  =  0),  Corresponding  to  equations 
(16)  different  expressions  and  obtained  for  the  positive  and  negative  ions 
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where  77  denot  the  location  of  the  peak  temperature  and  reaction 
'  max 

products  in  the  boundary  layer,  and  the  assumption  of  local  similarity 
eliminates  terms  on  the  right  hand  side  of  the  differential  equation.  In 
accordance  with  statement  (18)  pure  air  solutions  (e.g.,  references  1,2,  3, 
and  5)  can  be  used  to  determine  properties  of  function  F: 


10 


so  that  for  a  given  ionization  chemistry  (28)  represent  a  system  of  non¬ 
linear  algebraic  equations  for  the  unknown  peak  concentrations  of  the 

ionization  products,  c.  ,  in  the  boundary  layer.  The  thermodynamic 

^majc 

properties,  velocities,  concentration  of  major  (non-ionized)  species  are 

assumed  to  remain  unchanged  from  the  numerical  solution  obtained  for 

pure  air  injection.  The  total  contaminant  content  at  T7  is  known  as  a 

maoc 

function  of  injection  parameter  I  and  the  contaminants  concentration  'ji  the 
surface  material  To  complete  description  of  the  problem,  tfie  ionization 

chemistry  of  the  cesium-air  and  sodium-air  systems  must  now  be  described. 

The  system  of  chemical  reactions  that  was  considered  is  given  belcA^: 
Non-ionizihg  air  reactions: 

1)  0^*0^=  0*0*0^ 

2)  0  +0  =0  +  0  +  0 

3)  +  M  O  +  O  +  M  M  =  N,  N  ,  NO,  NO''’ 

4)  ^  N  +  N  +  N^ 

5)  N^  +  N  r  N  +  N+  N 

6)  N  +  M,  r?  N  +  N  +  M  =  O,  O^,  NO,  NO"^ 

D  6  6  2 

7)  NO  +  M  ^  N  +  O  +  M  M  =  O,  O  ,  N,  N  ,  NO^ 

I  I  I  w  M 

8)  O  +  N  +  NO 

9)  O  +  NO  N  +  O^ 

Pure  air  ionization 

10)  N  +  O  NO^  +  e’ 

11)  O  +  e’  O”  +  h 

V 

12)  O  +  e‘  +  -  o'  '  ^12  ^  ^2’ 

13)  02  +  e-  +  M,3.-0^-  +  Mj3  = 

The  air-cesium  and  air-sodium  ionization  reaction  systems  are 
identical  differing  only  in  their  rate  constants,  e.  g.  for  air  cesium  ionization 
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14)  Cs  +  O  =  Cs''"  +  o' 

W  M 

15)  Cs  +  O  ^  +  o' 

16)  Cs+M,,j=  Cs''"  +  e"  +  M,, 

io  lo  lo  Ic 

17)  Cs  +  e  n  Cs^  +  e  +  e 

18)  Cs^  +  e  s  Cs  +  h 

V 

19)  Cs  +  NO"^  ^  Cs''"  +  NO 

7 

The  rate  constants  for  reactions  1  through  10  are  given  in  Bortner's 
report;  they  are  identical  to  the  chemical  system  used  in  references  1,  Z 
and  5,  For  reactions  11  and  12,  and  the  sodium  system  (reactions  14 

g 

through  19)  Bortner's  latest  rate  constants  were  used.  Reaction  13  was 
assumed  to  be  in  equilibrium,  the  corresponding  equilibrium  constant,  and 
the  rate  constants  for  the  cesium  system  (reactions  14  through  19),  are 
based  on  Browne  and  Schwartz's^  work.  These  rate  constants  are  all 
listed  in  Appendix  I,  including  additional,  and  unrealistically  slow,  rate 
constants  for  reaction  12.  These  "slow"  rates  were  used  to  explore  the 
povential  significance  of  increased  O  concentrations  on  the  ionization 
process.  It  may  be  observed  that  the  five  unknown  charged  particle  con¬ 
centrations  can  be  reduced  to  a  system  of  three  unknowns  by  the  use  of 
overall  charge  conservation  and  the  partial  equilibrium  of  reaction  13. 

The  solution  of  the  ionization  problem  then  reduces  to  solving  non-linear 
algebiaic  system  (28)  for  three  charged  particle  concentrations,  arbitrarily 
chosen  as  Cs^  or  Na^,  NO^,  and  O  .  The  result  will  be  a  function  of 
distance  along  the  body,  the  thermo-chemical  state  of  the  uncontaminated 
air  at  tie  peaJc  temperature  point  in  the  boundary  layer,  and  the  total 
contamii  ant  content  at  that  point.  The  latter  depends  on  the  various  factors 
that  are  c'etermined  by  the  injection  process. 
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NUMERICAL  RESULTS 


Numerical  solutions  of  equations  {28)  were  obtained  for  conditions 
summarized  in  Table  I. 

TABLE  I 


Cone  half  angle 

PRINCIPAL  FLIGHT  CONDITIONS 

Altitude  U(n  Ug  Og 

k  ft,  kft, /sec  10"^  atm 

T 

^  w 
°K 

Contaminant 

10° 

150 

22 

21.  59 

27.56 

1000 

Cs 

10° 

200 

22 

21.  59 

4.6706 

1000 

Cs 

10° 

150 

22 

21.  59 

27.  56 

1000 

Na 

10° 

75 

20. 

72  20.334 

747.61 

1555.6 

Na 

9.03° 

150 

24. 

5  24.129 

27.616 

1450 

Na 

The  distances,  x,  from  the  tip  of  the  cone  ranged  from  0,  1  to  10  ft.  In 
addition  to  the  principal  flight  conditions  several  ’  hybrid"  conditions  were 
also  investigated  in  order  to  determine  the  effect  of  changing  only  one  of  the 
variables  (i,  e.  density,  temperature,  etc.)  on  the  solution.  It  was  found 
that  the  greatly  reduced  diffusivity  of  the  negative  ions  had  no  significant 
effect  in  reducing  peak  electron  densities,  except  for  the  more  slowly  ionizing 
sodium  case,  and  only  for  artificially  reduced  rates  of  the  O  detachment 
reaction  (reaction  12).  For  this  reason,  the  special  properties  of  the  negative 
ions  were  neglected  in  the  majority  of  the  calculations.  Table  II  presents 
a  summary  of  the  principal  results.  Since  the  foregoing  treatment  of  the 
ionization  problem  is  completely  separate  and  distinct  from  the  injection 
process.  Table  11  is  presented  in  terms  of  the  total  contaminant  content  at 

S  S' 

the  peak  temperature  ooint  (C^  (Tj  )  or  C  (77  )),  which  quantity 

max  max 

in  turn  depends  on  the  injection  parameters  (equations  23  and  26). 
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IONIZATION  AT  THE  PEAK  TEMPERATURE  POINT  IN  THE  BOUNDARY  LAYER 
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Examination  of  the  results  of  Table  II  indicated  that  only  the  forward 
ionization  rates  of  the  contaminant  were  significant  in  determining  the  fraction 
of  contam;.nant  that  was  ionised.  The  non-linear  system  (28)  can  then  be 
greatly  simplified  to  get  an  explicit  solution  for  the  fraction  of  sodium  or 
cesium  that  is  ionized  at  the  peak  temperature  point: 


c  max 
Na 

T 

c  max 


K  u 


1.0+ - - - 

^  /^14  ‘'O2 

- 32 -  ^ 


^15  ‘'O 


-1.  0 


Equation  (30)  was  found  to  represent  the  %  of  contaminant  ionized  (Table  III) 
within  5%.  Discrepancies  arose  only  at  the  highest  densities  of  charges  con¬ 
taminant  particles,  where  the  loss  due  reaction  18  amounted  to  a  few  percent. 
Reaction  15  was  found  to  be  the  most  rapid,  but  due  to  the  small  concentrations 
of  O,  reaction  (14)  often  dominated  the  result.  The  effect  of  O  concentration  is 
shown  explicitly  in  one  of  the  cases  in  Table  II  where  was  changed  from 
.  278  to  .  00278%.  It  is  expected  that  for  conditions  where  the  electron  densities 
are  significantly  larger  than  for  the  cases  presented  in  Table  II  reaction  17 
(ionization  of  contaminant  by  collision  with  electrons)  may  also  have  to  be  con¬ 
sidered.  In  that  case,  the  simple  formula  (30)  may  have  to  be  modified. 

Table  II  also  shows  the  close  relationship  between  the  particle  densities 

of  electrons  and  contaminant  ions;  differences  were  due  to  the  presence  of 

+  -  ^  ^ 

three  other  charged  particles;  NO  ,  O  and  .  For  all  cases  that  were 
considered,  the  presence  of  did  not  influence  electron  density  to  the  four 
significant  figures  that  are  presented  in  Table  II.  Only  for  the  cases  where 
the  total  contaminant  concentrations  were  very  small  did  the  presence  of  NO^ 
ions  add  significantly  to  the  peak  electron  density.  It  is  noteworthy  that  there 
was  no  effective  coupling  between  the  production  of  NO  ions  and  the  ionization 
of  the  contaminant.  The  reason  for  this  is  that  the  reactions  involving  charged- 
neutral  collisions,  which  may  have  coupled  the  two  types  of  ionization  processes. 
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occurred  much  less  frequently  than  the  ion  producing  neutral-neutral 
collisions.  The  only  detectable  change  in  NO  production  from  the  pure  air 
case  was  due  to  a  slight  depletion  of  the  NO^  producing  oxygen  atom  supply 
due  to  the  presence  of  O  . 

Equation^  (28)  and  (11)  imply  that  in  a  neutral  plasma  diffusion  of  the 
O  (and  other  negative)  ions  may  be  reduced  by  as  much  as  two  orders  of 
magnitude.  This  retardation  of  the  diffusion  is  due  to  the  microscopic  electric 
field  created  by  the  tendency  of  electrons  to  "escape”,  the  resulting  electrically 
attracting  region  then  hampers  diffusion  of  the  larger  negative  particles.  In 
the  ionizing  boundary  layer  this  retarded  diffusion  may  cause  an  "accumulation" 
of  O  (or  other  negative)  ions  in  the  reaction  zone,  greatly  reducing  thereby 
the  electron  densities  in  the  flow.  The  importance  of  this  effect  was  investi¬ 
gated  by  comparing  calculations  where  the  retarded  diffusivity  of  O  was  not 
accounted  for  with  the  more  exact  result.  An  unrealistically  slow  reaction 
rate  for  O  detachment  was  also  compared  to  the  exact  result.  The  results 

of  these  comparisons  are  shown  in  Table  III.  The  examples  for  this  com- 

o 

parison  were  taken  from  Table  II  for  sodium  contamination;  a  9:03  cone 
flying  at  24,  500  fps  and  a  10*^  cone  flying  at  22,  000  fps,  both  at  150,  000  ft. 
altitude  were  considered. 

TABLE  III 


The  results  of  Table  III  show  that  the  detachment  of  electrons  from  O  is  so 
rapid  that  the  slow  diffusion  rates  have  no  effect  on  the  outcome  of  the  calcu¬ 
lations.  At  the  "slow"  detachment  rates  the  outcome  is  quite  different,  the 
accumulation  of  O  significantly  reduced  peak  electron  densities  in  the 
boundary  layer.  This  indicates  that  it  may  be  important  for  some  types  of 
negative  ions  to  calculate  the  details  of  the  diffusion  process  in  reacting 
laminar  boundary  layers  in  order  to  predict  ionization  accurately.  This  is 
not  necessary  for  O  because  of  the  rapid  disappearance  of  this  ion  due  to 
electron  detachment  which  counteracts  the  ^  production  of  the  ionizing 
contaminant.  For  the  conditions  that  were  c'^nsidered,  reactions  15  (forward) 
and  12  (reverse)  were  the  only  two  reactions  that  affected  production  of  O  , 
so  that  it  was  possible  to  solve  equation  (28)  explicitly  for  the  peaPc  O  con¬ 
centration: 


■^15 


O 


max 


(k',  ^  ,  M)+  (.015-eK  u  /Sep  x) 
.2  e 


7?  =  T) 


max  (31) 


where  electron  and  positive  particle  concentrations  were  assumed  approxi' 
mately  equal,  so  tliat; 


& 


/n  + 

( /  lii 


m 

max  € 


^  C  max  k 


^  0.005 


(32) 


Equation  (30)  and  the  sodium  (cesium)  conservation  law  can  be  used  to 
rewrite  (31); 


O  max 


px 


c  )  c  ]  (G  +  px)  (G  +  ^ox) 

O  max  '  'max' 


(33) 


where ; 
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3X  K  u 


Sc 


"  ^14  ^O. 

_ _  ^ 

_  52 


+ 


k,^c 

15  o 

r6 


16 

“KT 


.  015  £  K  u 


Sc(k^^/23) 


12 


23 


ki^  M 


(34) 


Result  (33)  reaches  a  peak  value  when; 


P  X 


/ 


at  which  point; 


(35) 


O  max 


)  c^^(Tj  ) 
O  max  max 


2(GjG2G3) 


ITT 


(36) 


+  G  G  + 
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It  can  be  seen  that  when  G^  is  very  large  (fast  electron  detachment 
from  O  )  the  value  of  G^  is  immaterial,  because  the  O  concentration  will 
always  be  small,  approaching  in  the  limit.  Physically,  this  means 

that  in  this  limiting  case  the  concentration  of  O  is  determined  by  an 
"equilibrium"  between  the  production  of  O  by  reaction  15  forward  and  the 
consumption  of  O  by  reaction  12  in  reverse;  diffusion  plays  no  role.  For 
G^  not  large,  a  small  G^  (reduced  diffusivity)  will  increase  the  peak  con¬ 
centrations  of  the  negative  ions  greatly.  This  explains  the  results  of 
Table  III. 


18 


DISCUSSION,  CONCLUSIONS 


The  applicability  of  the  approximate  method  that  has  been  used  in  the 
foregoing  calculations  is  limited  by  the  underlying  assumptions.  These 
include  the  requirement  that  all  reaction  products  originate  in  the  high  tem¬ 
perature  region  of  the  boundary  layer,  that  the  wall  be  catalytic,  and  that  the 
boundary  conditions  change  slowly  if  at  all,  to  insure  local  similarity.  The 
latter  implies  that  if  the  body  is  not  a  slender  cone  it  must  be  a  pointed 
slender  body  with  slowly  varying  surface  slope.  The  accuracy  of  the  method 
is  further  limited  by  the  fact  that  due  to  the  differences  between  the  chemical 
and  transport  properties  of  the  different  species  the  maxima  in  temperature 
and  the  various  reaction  products  do  not  exactly  coincide,  and  that  the 
correlation  factor  K  (eq.  29)  varies  from  species  to  species.  For  NO^  and 
Blottner's  pure  air  cases  (ref.  1)  K  was  found  to  vary  between  3.  1  and  3.  5. 
(This  was  also  true  of  the  other,  iion-ionized  reaction  products  with  the 
exception  of  N.  )  For  the  cesium  calculations  (Table  II.  ),  K  was  taken  as 
3.43,  for  the  sodium  cases  (which  varied  ever  a  wider  range  of  conditions) 
a  value  of  K  =  3.  3  was  used,  I  being  .  63  and  Sc  =  .  5  in  both  cases.  For 
large  injection  rates,  the  thickening  of  the  boundary  layer  will  make  K 
smaller.  Examination  of  a  few  preliminary  numerical  boundary  layer 
solutions  (Blottner,  ref.  1)  with  injection  suggests  that 

K  -  ^  injection)  ^  3. 3 

V  T? 

max  max 

may  adequately  describe  this  effect.  Locally  similar  solutions  with  injection 
have  been  examined  resulting  in  a  relation: 

r?  =  1.0  +  V  x  (ft.)  (38) 

max 

that  may  be  useful.  Relations  (3?)  and  (38)  have  not  been  used  in  the  fore¬ 
going  results,  and  must  be  considered  conjectural.  They  may  however  form 
the  basis,  together  witn  equation  (30),  of  correlating  results  of  exact  numerical 
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calculations,  as  they  become  available.  It  is  expected  that  K  will  also  be  a 
weak  function  of  pressure,  wall  temperature,  maximum  temperature,  and 
other  parameters  of  the  problem.  In  particular,  as  the  ionization  approaches 
100%  and  equilibrium  is  approached  the  shape  of  the  profiles  will  certainly 
change  toward  a  K  decreased  severalfold,  and  the  peak  ion  contaminant  con¬ 
centration  will  shift  its  location  toward  the  wall.  The  present  method  is  then 
not  well  suited  to  conditions  near  ionization  equilibrium. 

The  treatment  of  the  contaminant  injection  in  the  foregoing  simplified 
analysis  is  entirely  separate  from  the  ionization  problem.  It  is  based  on 
constant  (or  slowly  varying)  injection  velocity  and  contaminant  content  in  the 
ablating  material.  If  these  assumptions  are  not  true,  but  the  total  con¬ 
taminant  content  at  the  peak  temperatur ’■  point  can  be  calculated  or  estimated 
by  some  other  method,  the  calculation  of  the  ionization  can  still  proceed  by 
means  of  equations  (28).  For  most  of  the  results  that  were  calculated 

—  (Sc  T)  )  =  0.24 
u  max 

e 

=  0.40  (30) 

T)  =  0 

represent  reasonable  approximations. 

Finally,  comparison  of  the  approximate  method  uith  a  preliminarv 
exact  result,  heretofore  unpublished,  v^as  made  possible  thanks  to  the 
courtesy  of  the  author's  colleague  Dr.  F.  G.  Blottner.  Injection  of  sodium 
contaminated  air  was  consid^^  red  for  a  simplified  chemical  system  consisting 
of  reaction  1  through  12  and  15.  (The  approximate  results  were  als  com¬ 
puted  for  this  simplified  chemical  system;  i.e.  reactions  10  through  12  and  15.) 
The  conditions  of  the  comparison  and  the  re..ults  are  summarized  in  Table  IV. 


1  5  U 
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TABLE  IV. 


COMPARISON  OF  EXACT  AND  APPROXIMATE  SOLUTIONS 


Cone  angle 
Altitude 


T 

T 

max 

^(Tjmax) 


Input 


10° 

150,  000  ft 
22. 000  fps 
21,  5^0  fps 
1, 000°K 
4.478°K 

.2161  X  10  g/cc 
. 62845 

3.  1273 

4.  5809 

.  206  ^ 

.  33555  X  10 

1.  0  X  10‘^ 

1.2  ft 


Results  at  Tf 


Exact 


^Na  ,^5 

C  X  10 


10 

narticle  s 


1.  362 


45.  70 


.  5165 


0976 


46.  12 


Approximate 


1.  361 


40.  28 


.4164 


.  1097 


40.  59 


The  agreement  for  C^  up  to  the  fourth  significant  figure  must  be  considered 

fortuitous.  The  remaining  results  however  agree  within  the  expected  accuracy 

of  the  thermochemical  properties  of  the  reacting  gas. 

9 

At  the  conclusion  of  this  work  Kane's  report  for  sodium  contaminated 
air  ionization  came  to  the  author's  attention.  In  comparison  with  the  present 

9  .  .  ■  • 

work,  Kane  makes  the  same  series  of  fluid  mechanic  as  sumptions.  His 

treatment  of  the  injection  is  identical  except  for  a  more  complicated  form  of 

the  "Reynolds  analogy".  The  solution  of  the  equation  is  performed  by 

9 

numerically  evaluating  integrals  over  the  entire  boundary  layer  (Kane  did 
not  take  advantage  of  the  short  cut  afforded  by  equations  (28)  of  the  present 
report).  Finally,  Kane's  kinetics  for  sodium  are  taken  from  Bortner's 

earlier  work,  which  has  been  superseded  since.  (The  complicating  effects 

-  -  9 

of  O  ,  O  are  not  accounted  for  by  Kane  ,  this  may  not  result  in  serious 

2  9 

error.)  In  particular,  the  reaction  that  dominates  Kane's  results 


Na  +0+0 


Na  +  e  +  O, 
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should  be  replaced  by  reaction  14  and  15.  Finally,  Kane  does  not  account 

for  the  effect  of  d.epletion  of  the  contaminant  material  on  the  ionization  rate. 

9  + 

Thus,  Kane's  equation  (36)  gives  a  linear  growth  law  for  N  vs  x,  which 

IN  ^ 

is  equivalent  to  expanding  equation  (30)  of  the  present  analysis  in  powers  of 
X.  The  more  accurate  expression,  eq.  (30),  is  superior  when  the  fraction 
of  contaminant  ioni.’.ed  becomes  significant. 

In  conclusioi  ,  it  can  be  stated  that,  subject  to  the  assumptions  of 
the  analysis: 

(1)  The  apprc'ximate  method  gives  reasonable  predictions  of  the 
non-equilibrium  ionization  process  in  contaminated  laminar  hypersonic 
boundary  layers  on  slender  bodies. 

(2)  The  injection  and  ionization  processes  are  independent. 

(3)  The  ionizations  of  pure  air  and  the  alkali  metal  contaminant 
are  decoupled. 

(4)  The  ionization  of  alkali  metal  contaminants  in  laminar  hypersonic 
slender  body  boundary  ’  yers  is  dominated  by  the  forward  ionization  rates 
and  can  be  expressed  by  a  simple  formula  (eq.  30),  unless  contaminant 
ionization  due  to  electron -contaminant  collision  (reaction  7)  becomes 
important. 

(5)  Due  to  the  high  detachment  rate  of  electrons  from  O  ,  the 
presence  of  these  ions  does  not  affect  the  result  significantly  in  spite  of  their 
retarded  diffusion  rates. 

(6)  The  results  of  this  approximate  method  can  be  used  to  make  pre¬ 
dictions  of  contaminant  ionization  and  to  correlate  and  explain  results  of  exact 
calculations  as  they  become  available. 
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NOMENCLATURE 


c. 

1 


D.. 

ij 


Mass  fraction  of  i-th  species 

Mass  fraction  of  i-th  element 

Mixture  diffusion  coefficient 

Binary  diffusion  coefficient 

Body  force  due  to  microscopic  electric  field 


f  Nor*-dimensional  stream  function  (Equation  2) 

F  Similarity  function  (Equation  18) 

Equation  (34) 

h  Mixture  enthalpy  • 


h. 


K 


m 


M 


n 


n  n 
P  n 

N 


Species  enthalpy 

Injection  parameter  (Equation  27) 

Mixture  (frozen)  heat  conductivity 

Chemical  rate  constants  for  11-th,  12-th  etc.  reaction 


K  =  -a  F/ 


(Equation  29) 


max 


Density-viscosity  product  (Equation  5) 
Molecular  weight  of  i-th  species 
Molecular  weight  of  mixture 
Number  of  species  in  mixture 


Number  of  positive  (negativ*)  ion  species  in  mixture 

Number  densities 

Pressure 


Pr 


(Frozen)  Prandtl  number  (Equation  5) 


24 


Diffusional  flux  velocity  of  i-th  species 
r  Radial  coordinate  of  axi-symmetric  surface 

R  Universal  gas  constant 

Sc  Schmidt  number 

T  Absolute  temperature 

u,  V  Parallel  and  normal  velocities  in  boundary  layer 

w.  Chemical  mass  generation  rate  of  i-th  species 

X,  y  Parallel  and  normal  coordinates  in  boundary  layer 

Number  of  atoms  of  element  j  in  species  i 
j8  Pressure  gradient  parameter  (Equation  5) 

^  I  T]  Transformed  parallel  and  normal  coordinates  in  boundary  layer 

(Equation  1) 

M  Viscosity 

p  Density 

SUBSCRIPTS 

e  At  the  outer  edge  of  the  boundary  layer 

f  Electrons 

eq  In  chemical  equilibrium 

i,j,0  Refers  to  species  (i-th,  j-th,  O,  etc.) 

max  Maximum,  or  at  maximum  temperature 

w  W  all 

”  Free  stream 

SUPERSCRIPTS 

i,j,0  Refers  to  element  (i-th,  j-th,  O,  etc.)  _ 

,  -  Charged  particles 
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APPENDIX 

CHEMICAL  RATE  CONSTANTS 

The  units  are  °K  (temperature),  g  moles /cc  (concentration)  and 

cm  /gmole 

b  c 

sec  or  cm  /(  g  mole)  sec  (rate  constants  for  2  and  3  body 

reactions  respectively).  Also;  R  =  1,9864  cal/°Kmole 

'  • 

Reaction 

"^FORWARD 

^REVERSE 

X 

1.9 

xlo'S-'/^ 

2 

7.  1 

x,0‘S-‘'' 

3 

,  n  -59.400/T 

3.0x10  T  e 

2.  5 

XIO'^T-*/^ 

4 

2  X 

10’«T-‘ 

5 

7.x 

6 

1.9xl0'\-'e-“5-2»»/T 

1  X 

,  .  1  ‘3  -  1 

10  T 

7 

.  ,  ,.17  -1  /2  -75,  500/T 

2.4x10  T  e 

6  X 

,.16  -1/2 

10  T 

8 

,  O  ,nl3  -37,  750/T 

6.  8  X  1 0  e 

1.5 

X  10^^ 

9 

7  3/2  -19, 100/T 

4.3x10  T  e 

1.8 

10 

1  ^  -31,900/T 

1.3x10  T  e 

2  X 

1  1 

7.2  X  10® 

•  . 

12 

19  -  5 

8. 5  X  10  T 

1..  4 

,^12^  -33800/RT 

X  10  T  e 

1  3 

,  _  ^  ,.-9_1.44  -10650/RT 

k  x7,  3xl0T  e 

reverse 

Cesium 

14 

,^17^.07  -79130/RT 

1 0  T  e 

2  X 

10^^ 

15 

_  _  ,.16  15  -55970/RT 

8.7x10  T  e 

2  X 

io^'^t"*  ® 

16 

,  -  ,.12„1.  18  -89730/RT 

1,2x10  T  e 

3  X 

lo'^ 

17 

,  ,  ,.32^-3.  3  -89780/RT 

1.6x10  T  e 

4  X 

18 

_  _  ,.14  -.67 

7.2x10  T 

• 

19 

13  -  5 

6.  02  X  10  T 

807  xlo‘‘e-‘22^"“«T 

Sodium 

14 

,  „  ,„15  -108260/RT 

1.8x10  e 

1.2 

15 

^  ,^15  -854]5/RT 

2.0x10  e 

1.2 

xl0'\-' 
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Reaction 

FORWARD 

REVERSE 

16 

3.9 

,^6  1.  5  -1  18590/RT 

X  1  0  T  e 

9.  8 

14 

X  10 

17 

3.9 

,„14  .  5  -11S590/RT 

X  1 0  T  e 

00 

• 

CT' 

22  -1.0 

X  IG  T 

18 

i.2 

19 

3  X 

9  X 

,^10  -94753/RT 
10  e 

12  slow 

2.  7 

5  X 

,-5^2  -33800/RT 
1 0  T  e 
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SUMMARY 

An  approximate  procedure  is  used  to  predict  the 
effect  on  air  ionization  of  small  amounts  of  cesium  and 
sodium  in  the  ablating  surface  material  of  pointed  cones. 
Ionization  is  assumed  to  occur  due  to  finite  rate  gas  phase 
chemical  reactions  in  the  laminar  boundary  layer.  An 
11  species  (N^#  0^$  N,  O,  NO,  Cs  or  Na,  NO"^,  Cs+  or 
Na'*',  0_',  O",  e")  16  reaction  chemical  system  is 
assumed. 
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